33411  2 


• l*o 

reports  Control  Symbol 
OSD-1366 


F.COM-5W6 


THE  DEPTH  OF  THE  SURFACE 
BOUNDARY  LAYER 


By 

Frank  V.  Hansen 


Atmospheric  Scisncss  Laboratory 


Juno  1976 


UNITEO  STATES  ARMY  ELECTRONICS  COMMAND  - FORT  MONMOUTH  TW  JERSEY  07703 


NATIONAL  TEOMCAL 
^FORMATION  SERVICE 


u *■  SSmimS.  * aw 


ft 


• V 


NOTICES 


Disclaimers 

The  findings  in  this  report  are  not  to  be  construed  as  an  of- 
ficial Department  of  the  Army  position,  unless  so  designated 
by  other  authorized  documents. 


The  citation  of  trade  names  and  names  of  manufacturers  in 
this  report  is  not  to  be  construe d as  official  Government  in- 
dorsement or  approval  of  commercial  products  or  services 
referenced  herein. 


Disposition 

Destory  this  report  when  it  is  no  longer  needed.  Do  not 
return  it  to  the  originator. 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (Whan  Data  Entarad) 


REPORT  DOCUMENTATION  PAGE 


. REPORT  NUMBER 

ECOM-5596 


4.  TITLE  (and  Subtltla) 


12.  GOVT  ACCESSION  NO. 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


3.  RECIPIENT’S  CATALOG  NUMBER 


s.  type  of  report  4 period  covereo 


THE  DEPTH  OF  THE  SURFACE  BOUNDARY  LAYER 


6.  PERFORMING  ORG.  REPORT  NUMBER 


7.  AUTHORS 


8.  CONTRACT  OR  GRANT  NUMBERO) 


Frank  V.  Hansen 


9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Atmospheric  Sciences  Laboratory 

White  Sands  Missile  Range,  New  Mexico  88002 


It.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

US  Army  Electronics  Command 
Fort  Monmouth,  New  Jersey  07703 


i 


MONITORING  AGENCY  NAME  8 ADDRESSfif  dIHarant  from  Controlling  Oltlca) 


10.  PROGRAM  ELEMENT.  PROJECT,  TASK 
AREA  ft  WORK  UNIT  NUMBERS 


DA  Task  No.  1T162111AH71A10 


t2.  REPORT  DATE 

June  1976 


13  NUMBER  OF  PAGES 

IS 


15.  SECURITY  CLASS,  (ot  thla  raport) 


UNCLASSIFIED 


IS*.  OECL  ASSI FIC  ATION/  DOWNGRADING 
SCHEDULE 


■ft.  DISTRIBUTION  STATEMENT  (ot  thla  Raport) 


Approved  for  public  release;  distribution  unlimited. 


17.  DISTRIBUTION  STATEMENT  (ot  (h«  «bi(r«c(  entered  in  Block  20,  If  different  from  Report ) 


19.  KEY  WORDS  (Continue  on  reeeree  e/de  If  neceeeery  end  Identify  by  block  number) 


Surface  boundary  layer 

Micrometeorology 

Dynamic  similarity  theory 


Richardson  number 
Reynolds  stress 


20.  ABSTRACT  (Caatmua  an  rararaa  aMa  M nmaaaaaty  mod.  Identity  by  block  nuatbar) 

The  depth  of  the  surface  boundary  layer  of  the  atmosphere  can  be  related  to  a 
characteristic  dimension  defined  by  Obukhov  as  the  scaling  length  L.  Utilizing 
the  form  of  the  diabatic  wind  and  temperature  profiles,  the  depth  of  the  surface 
layer  is  demonstrated  to  be  proportional  to  and  a function  of  the  slope  of  the 
wind  profile  and  the  Obukhov  length.  The  methodology  was  verified  utilising 
experimental  data  extracted  from  the  literature.  Generally,  the  surface  bound- 
ary layer  depth  approaches  the  depth  of  the  planetary  boundary  layer  in  near 
adiabatic  conditions  near  sunrise  and  sunset. 


FORM 
I JAM  73 


EDITION  OF  • NOV  «S  IS  OBSOLETE 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  fWB«n  Data  Entarad) 


SUmARY 


In  general,  the  approach  presented  to  determine  the  depth  cf  the  sur- 
face boundary  layer  Is  simple  and  relatively  straightforward.  The  depths 
of  the  surface  boundary  layer  as  indicated  in  Figure  1 are  representative 
averages,  while  the  depths  shown  in  the  remaining  figures  are  relative 
owing  to  the  smoothing  techniques  utilized.  The  binomial  filtering 
method  has  a tendency  to  reduce  peaks  drastically.  The  results  obtained 
appear  to  be  reasonable  and  agree  with  the  experimental  determinations 
of  Lumley  and  Panofsky  . 


CONTENTS 

Page 


SUWARY  1 
INTRODUCTION  3 
DYNAMIC  SIMILARITY  IN  THE  SURFACE  BOUNDARY  LAYER  3 
THE  DEPTH  OF  THE  SURFACE  BOUNOARY  LAYER  5 
DISCUSSION  7 
REFERENCES  8 


2 


or* . "j*.  wi 

a 


T 


'£ 


w, 


INTRODUCTION 


Within  the  planetary  boundary  layer  known  as  the  friction  zone  of  the 
atmosphere,  a sublayer  exists  where  the  heat  and  momentum  fluxes  are 
considered  to  be  Invariant  with  height.  This  surface  boundary  layer 
is  further  characterized  by  a scalar  Reynolds  stress  along  the  mean 
wind,  the  direction  of  which  Is  also  invariant  with  height  since  Cori- 
olis effects  may  be  assumed  to  be  Insignificant.  The  depth  of  the  sur- 
face boundary  layer  can  vary  from  a few  meters  to  perhaps  a hundred  or 
more  meters,  depending  upon  surface  roughness,  stability,  and  windspeed. 
In  turn,  these  parameters  and  the  thickness  of  the  surface  boundary  layer 
control  the  maximum  values  of  the  eddy  conductivity  and  eddy  viscosity 
that  occur  in  the  planetary  boundary  layer. 


A method  for  determining  the  thickness  (or  depth)  of  the  surface  boundary 
layer  can  be  easily  developed  from  the  dynamic  similarity  theory  of 
Obukhov  [1],  in  terms  of  measurable  and  calculated  micrometeorological 
parameters.  The  scheme  is  simple,  readily  implemented,  and  cllows  rapid 
determination  of  the  depth  of  the  surface  boundary  layer  and  exchange 
coefficient  maxims  for  the  planetary  boundary  layer. 


DYNAMIC  SIMILARITY  IN  THE  SURFACE  BOUNDARY  LAYER 


Obukhov  [1]  postulated  that  heat  and  momentum  fluxes  in  the  surface 
boundary  layer  were  rynamically  similar.  Drawing  on  Richardson's  [2] 
criterion  and  the  Schmidt  [3]  exchange  coefficient  hypothesis,  Obukhov 
reasoned  that  heat  and  momentum  transfer  processes  in  the  surface  boundary 
layer  could  be  stated  as 
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where  H is  the  vertical  heat  flux,  c the  specific  heat  of  air  at  con- 
stant pressure,  e,  the  potential  temperature,  V the  mean  windspeed, 
the  ratio  of  the  exchange  coefficients  for  heat  and  momentum, 

p density,  x the  Reynolds  stress,  0 (Ri)  a monotonously  decreasing  uni- 
versal function,  i the  Prandtl  [4]  mixing  length,  an-,  z height.  Paral- 
leling, and  In  conjunction  with  Equations  (1)  and  (2),  the  conjugate 
laws  for  the  wind  and  temperature  profiles  In  differential  form  may 
be  written  as: 
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where  u*  is  a friction  velocity,  k Karman's  constant,  0M  = [0(Ri)]"\ 
0H  = [0(Ri)]_1,  and  T*  a scaling  temperature  defined  as 


ku*  cpp 
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The  gradient  Richardson  number,  Ri , is  generally  written  as 


Ri  = a je/». 
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where  g is  the  gravitational  acceleration.  Operating  upon  Equation  (8) 
with  respect  to  0(Ri),  Obukhov  found  that  Equation  (8)  could  be  stated 
as 
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If  Equation  (9)  is  differentiated  with  respect  to  z,  what  is  now  known 
as  the  Obukhov  scaling  length  is  represented  by 
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where  L is  a characteristic  length  representative  of  and  proportional 
to  the  depth  of  the  surface  boundary  layer.  By  definition  it  may  also  be 
stated  that 
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which  may  also  be  stated  as  z/L  = (L  Rf  where  R.  is  the  flux  form  of  the 
Richardson  number. 
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A universal  form  of  the  wind  profile  nuy  now  be  developed  from  Equations 
(3)  and  (11).  Writing  Equation  (3)  as 
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Then  adding  and  subtracting  1 to  z/L  R*  1 and  multiplying  and  dividing 
by  z/L, 


Defining  z/L  - R.(Rx 
Equation  (13),  T T 
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z/L)"1  as  an  arbitrary  variable  6,  then  integrating 
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where  6 is  the  average  8 over  the  layer  zg 
roughness  length.  Since  Rf  = Ri  l^/K^,  then 
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which  is  valid  for  both  stable  and  unstable  thermally  stratified  flow 
in  the  surface  boundary  layer. 


THE  DEPTH  OF  THE  SURFACE  BOUNDARY  LAYER 


Lumley  and  Panofsky  [5]  suggest  that  the  depth  of  the  surface  boundary 
layer  may  be  enplrfcally  determined  from 


where  h is  In  meters  and  tq  is  the  surface  stress.  It  may  be  assumed 
that  t0  3 t * p u£,  the  Reynolds  stress.  Equation  (16)  is  based  upon 
the  assumption  that  t varies  about  20  percent  from  zQ  to  h. 

A more  precise  argument  can  be  developed  from  Equation  (3).  If  Equation 
(3)  is  written  in  finite  difference  form  fo*  the  unstable  regime  as 


where  z = h = L Ri , then 


If  aV  the  wind  gradient  at  the  geometric  mean  height  h is  assumed  to 
occur  over  the  layer  Ain  z = In  e = 1,  then  AV  = Ri  and 


the  slope  of  the  normalized  diabatic  profile. 
In  thermally  stratified  stable  flow.  Equation  (3)  is  written  as 


where  m 


owing  to  the  assumption  that  all  fluxes  are  a result  of  mechanical 
turbulence  only.  This  leads  to  and  z/L  = Ri  0^.  If  Ain  z is 

again  unity,  then  Ri  aV  = 8-1  and 


since  for  stable  flow,  Hansen  [6]  has  shown  that 


If  h/L  - R1  0^,  then  from  Equation  (21 ) 


R*i  = Uj 

Rl(h)  15 


* L[R1(h) + 15 


whore  the  subscriot  denotes  the  height  of  the  parameter  evaluation. 

DISCUSSION 

Equations  (19)  and  (21),  the  primary  formulae  for  calculating  the  depth 
of  the  surface  boundary  layer,  were  evaluated  by  using  experimental  data 
extracted  from  studies  reported  on  by  Lettau  and  Davidson  [7],  Barad  [8], 
Swinbank  [9]  and  Izumi  [10].  Generally,  these  data  were  observed  in 
relatively  stationary  conditions  and  mostly  over  terrain  that  possessed 
a high  degree  of  homogeneity.  Of  the  493  profiles  available,  only  36  were 
unusable  because  of  the  extremely  low  windspeeds  reported. 

Gradient  Richardson  numbers,  scaling  lengths,  profile  slopes,  Reynolds 
stresses,  and  the  dimensionless  shears  were  calculated  for  each  profile. 
Surface  boundary  layer  depths  were  determined  for  both  the  stable  and 
unstable  flow  regimes.  These  data  were  then  averaged  as  a geometric 
progression  as  a function  of  L and  are  shown  in  Figure  1.  The  curve 
representing  h,  calculated  by  using  Equation  (16),  is  also  given.  The 
agreement  between  Equation  (16)  and  h as  determined  from  Equations  (19) 
and  (21)  is  considered  to  be  good.  Included  for  definitive  purposes  are 
the  smoothed  average  windspeeds  observed  over  the  entire  stability  range. 

The  divergence  of  the  two  curves  over  the  range  -55  < L < 222  meters  is 
attributed  to  the  fact  that  the  profiles  observed  in  this  stability  range 
were  observed  near  sunrise  and  sunset  under  light  wind  conditions  and 
indifferent  thermal  stratification.  The  dashed  portions  of  the  curves 
through  neutral  conditions  are  suggested  shapes  owing  to  a lack  of  data 
near  L = infinity. 

Surface  boundary  layer  depth  with  respect  to  windspeed  at  2 m above  the 
surface  is  shown  in  Figure  2.  Both  V and  h were  smoothed  with  a five-point 
binomial  filter  over  24  hours  to  obtain  this  representation.  It  is  apparent 
that  surface  boundary  layer  depth  is  somewhat  of  a function  of  windspeed 
and.  In  turn,  the  Reynolds  stresses. 

Figures  3a  and  3b  present  the  surface  boundary  layer  depth  as  a function 

of  L and  L"1 , respectively.  Again,  these  data  were  smoothed  by  using  the 
binomial  filter  over  24  hours.  Figures  3a  and  3b  show  that  in  near  neutral 
conditions,  the  depth  of  the  surface  boundary  layer  approaches  the  depth 
of  the  planetary  boundary  layer.  If  this  Is  so,  then  the  limit  of  the 
Obukhov  [1]  length  will  not  be  Infinity  but  the  overall  depth  of  the  plane- 
tary boundary  layer  in  adiabatic  flow  conditions. 
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